A real-world data validation of the value of early-stage SIR modelling to public health.

Performance of Susceptible-Infected-Recovered (SIR) model in the early stage of a novel epidemic may be hindered by data availability. Additionally, the traditional SIR model may oversimplify the disease progress, and knowledge about the virus and transmission is limited early in the epidemic, resulting in a greater uncertainty of such modelling. We aimed to investigate the impact of model inputs on the early-stage SIR projection using COVID-19 as an illustration to evaluate the application of early infection models. We constructed a modified SIR model using discrete-time Markov chain to simulate daily epidemic dynamics and estimate the number of beds needed in Wuhan in the early stage of COVID-19 epidemic. We compared eight scenarios of SIR projection to the real-world data (RWD) and used root mean square error (RMSE) to assess model performance. According to the National Health Commission, the number of beds occupied in isolation wards and ICUs due to COVID-19 in Wuhan peaked at 37,746. In our model, as the epidemic developed, we observed an increasing daily new case rate, and decreasing daily removal rate and ICU rate. This change in rates contributed to the growth in the needs of bed in both isolation wards and ICUs. Assuming a 50% diagnosis rate and 70% public health efficacy, the model based on parameters estimated using data from the day reaching 3200 to the day reaching 6400 cases returned a lowest RMSE. This model predicted 22,613 beds needed in isolation ward and ICU as on the day of RWD peak. Very early SIR model predictions based on early cumulative case data initially underestimated the number of beds needed, but the RMSEs tended to decline as more updated data were used. Very-early-stage SIR model, although simple but convenient and relatively accurate, is a useful tool to provide decisive information for the public health system and predict the trend of an epidemic of novel infectious disease in the very early stage, thus, avoiding the issue of delay-decision and extra deaths.

Modeling the spread of COVID-19 in spatio-temporal context.

This study aims to use data provided by the Virginia Department of Public Health to illustrate the changes in trends of the total cases in COVID-19 since they were first recorded in the state. Each of the 93 counties in the state has its COVID-19 dashboard to help inform decision makers and the public of spatial and temporal counts of total cases. Our analysis shows the differences in the relative spread between the counties and compares the evolution in time using Bayesian conditional autoregressive framework. The models are built under the Markov Chain Monte Carlo method and Moran spatial correlations. In addition, Moran's time series modeling techniques were applied to understand the incidence rates. The findings discussed may serve as a template for other studies of similar nature.

A coupled awareness-epidemic model on a multi-layer time-varying network.

Social interactions have become more complicated and changeable under the influence of information technology revolution. We, thereby, propose a multi-layer activity-driven network with attractiveness considering the heterogeneity of activated individual edge numbers, which aims to explore the role of heterogeneous behaviors in the time-varying network. Specifically, three types of individual behaviors are introduced: (i) self-quarantine of infected individuals, (ii) safe social distancing between infected and susceptible individuals, and (iii) information spreading of aware individuals. Epidemic threshold is theoretically derived in terms of the microscopic Markov chain approach and the mean-field approach. The results demonstrate that performing self-quarantine and maintaining safe social distance can effectively raise the epidemic threshold and suppress the spread of diseases. Interestingly, individuals' activity and individuals' attractiveness have an equivalent effect on epidemic threshold under the same condition. In addition, a similar result can be obtained regardless of the activated individual edge numbers. The epidemic outbreak earlier in a situation of the stronger heterogeneity of activated individual edge numbers.

An AI approach for managing financial systemic risk via bank bailouts by taxpayers.

Bank bailouts are controversial governmental decisions, putting taxpayers' money at risk to avoid a domino effect through the network of claims between financial institutions. Yet very few studies address quantitatively the convenience of government investments in failing banks from the taxpayers' standpoint. We propose a dynamic financial network framework incorporating bailout decisions as a Markov Decision Process and an artificial intelligence technique that learns the optimal bailout actions to minimise the expected taxpayers' losses. Considering the European global systemically important institutions, we find that bailout decisions become optimal only if the taxpayers' stakes exceed some critical level, endogenously determined by all financial network's characteristics. The convenience to intervene increases with the network's distress, taxpayers' stakes, bank bilateral credit exposures and crisis duration. Moreover, the government should optimally keep bailing-out banks that received previous investments, creating moral hazard for rescued banks that could increase their risk-taking, reckoning on government intervention.

New Surveillance Metrics for Alerting Community-Acquired Outbreaks of Emerging SARS-CoV-2 Variants Using Imported Case Data: Bayesian Markov Chain Monte Carlo Approach.

BACKGROUND: Global transmission from imported cases to domestic cluster infections is often the origin of local community-acquired outbreaks when facing emerging SARS-CoV-2 variants. OBJECTIVE: We aimed to develop new surveillance metrics for alerting emerging community-acquired outbreaks arising from new strains by monitoring the risk of small domestic cluster infections originating from few imported cases of emerging variants. METHODS: We used Taiwanese COVID-19 weekly data on imported cases, domestic cluster infections, and community-acquired outbreaks. The study period included the D614G strain in February 2020, the Alpha and Delta variants of concern (VOCs) in 2021, and the Omicron BA.1 and BA.2 VOCs in April 2022. The number of cases arising from domestic cluster infection caused by imported cases (Dci/Imc) per week was used as the SARS-CoV-2 strain-dependent surveillance metric for alerting local community-acquired outbreaks. Its upper 95% credible interval was used as the alert threshold for guiding the rapid preparedness of containment measures, including nonpharmaceutical interventions (NPIs), testing, and vaccination. The 2 metrics were estimated by using the Bayesian Monte Carlo Markov Chain method underpinning the directed acyclic graphic diagram constructed by the extra-Poisson (random-effect) regression model. The proposed model was also used to assess the most likely week lag of imported cases prior to the current week of domestic cluster infections. RESULTS: A 1-week lag of imported cases prior to the current week of domestic cluster infections was considered optimal. Both metrics of Dci/Imc and the alert threshold varied with SARS-CoV-2 variants and available containment measures. The estimates were 9.54% and 12.59%, respectively, for D614G and increased to 14.14% and 25.10%, respectively, for the Alpha VOC when only NPIs and testing were available. The corresponding figures were 10.01% and 13.32% for the Delta VOC, but reduced to 4.29% and 5.19% for the Omicron VOC when NPIs, testing, and vaccination were available. The rapid preparedness of containment measures guided by the estimated metrics accounted for the lack of community-acquired outbreaks during the D614G period, the early Alpha VOC period, the Delta VOC period, and the Omicron VOC period between BA.1 and BA.2. In contrast, community-acquired outbreaks of the Alpha VOC in mid-May 2021, Omicron BA.1 VOC in January 2022, and Omicron BA.2 VOC from April 2022 onwards, were indicative of the failure to prepare containment measures guided by the alert threshold. CONCLUSIONS: We developed new surveillance metrics for estimating the risk of domestic cluster infections with increasing imported cases and its alert threshold for community-acquired infections varying with emerging SARS-CoV-2 strains and the availability of containment measures. The use of new surveillance metrics is important in the rapid preparedness of containment measures for averting large-scale community-acquired outbreaks arising from emerging imported SARS-CoV-2 variants.

Stochastic Lag Time Parameterization for Markov State Models of Protein Dynamics.

Markov state models (MSMs) play a key role in studying protein conformational dynamics. A sliding count window with a fixed lag time is widely used to sample sub-trajectories for transition counting and MSM construction. However, sub-trajectories sampled with a fixed lag time may not perform well under different selections of lag time, which requires strong prior practice and leads to less robust estimation. To alleviate it, we propose a novel stochastic method from a Poisson process to generate perturbative lag time for sub-trajectory sampling and utilize it to construct a Markov chain. Comprehensive evaluations on the double-well system, WW domain, BPTI, and RBD-ACE2 complex of SARS-CoV-2 reveal that our algorithm significantly increases the robustness and power of a constructed MSM without disturbing the Markovian properties. Furthermore, the superiority of our algorithm is amplified for slow dynamic modes in complex biological processes.

Novel Type I Half Logistic Burr-Weibull Distribution: Application to COVID-19 Data.

In this work, we presented the type I half logistic Burr-Weibull distribution, which is a unique continuous distribution. It offers several superior benefits in fitting various sorts of data. Estimates of the model parameters based on classical and nonclassical approaches are offered. Also, the Bayesian estimates of the model parameters were examined. The Bayesian estimate method employs the Monte Carlo Markov chain approach for the posterior function since the posterior function came from an uncertain distribution. The use of Monte Carlo simulation is to assess the parameters. We established the superiority of the proposed distribution by utilising real COVID-19 data from varied countries such as Saudi Arabia and Italy to highlight the relevance and flexibility of the provided technique. We proved our superiority using both real data.

Scalable Bayesian phylogenetics.

Recent advances in Bayesian phylogenetics offer substantial computational savings to accommodate increased genomic sampling that challenges traditional inference methods. In this review, we begin with a brief summary of the Bayesian phylogenetic framework, and then conceptualize a variety of methods to improve posterior approximations via Markov chain Monte Carlo (MCMC) sampling. Specifically, we discuss methods to improve the speed of likelihood calculations, reduce MCMC burn-in, and generate better MCMC proposals. We apply several of these techniques to study the evolution of HIV virulence along a 1536-tip phylogeny and estimate the internal node heights of a 1000-tip SARS-CoV-2 phylogenetic tree in order to illustrate the speed-up of such analyses using current state-of-the-art approaches. We conclude our review with a discussion of promising alternatives to MCMC that approximate the phylogenetic posterior. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Evaluating Effects of Dynamic Interventions to Control COVID-19 Pandemic: A Case Study of Guangdong, China.

The emergence of different virus variants, the rapidly changing epidemic, and demands for economic recovery all require continual adjustment and optimization of COVID-19 intervention policies. For the purpose, it is both important and necessary to evaluate the effectiveness of different policies already in-place, which is the basis for optimization. Although some scholars have used epidemiological models, such as susceptible-exposed-infected-removed (SEIR), to perform evaluation, they might be inaccurate because those models often ignore the time-varying nature of transmission rate. This study proposes a new scheme to evaluate the efficiency of dynamic COVID-19 interventions using a new model named as iLSEIR-DRAM. First, we improved the traditional LSEIR model by adopting a five-parameter logistic function ß(t) to depict the key parameter of transmission rate. Then, we estimated the parameters by using an adaptive Markov Chain Monte Carlo (MCMC) algorithm, which combines delayed rejection and adaptive metropolis samplers (DRAM). Finally, we developed a new quantitative indicator to evaluate the efficiency of COVID-19 interventions, which is based on parameters in ß(t) and considers both the decreasing degree of the transmission rate and the emerging time of the epidemic inflection point. This scheme was applied to seven cities in Guangdong Province. We found that the iLSEIR-DRAM model can retrace the COVID-19 transmission quite well, with the simulation accuracy being over 95% in all cities. The proposed indicator succeeds in evaluating the historical intervention efficiency and makes the efficiency comparable among different cities. The comparison results showed that the intervention policies implemented in Guangzhou is the most efficient, which is consistent with public awareness. The proposed scheme for efficiency evaluation in this study is easy to implement and may promote precise prevention and control of the COVID-19 epidemic.

Effects of the information-driven awareness on epidemic spreading on multiplex networks.

In this study, we examine the impact of information-driven awareness on the spread of an epidemic from the perspective of resource allocation by comprehensively considering a series of realistic scenarios. A coupled awareness-resource-epidemic model on top of multiplex networks is proposed, and a Microscopic Markov Chain Approach is adopted to study the complex interplay among the processes. Through theoretical analysis, the infection density of the epidemic is predicted precisely, and an approximate epidemic threshold is derived. Combining both numerical calculations and extensive Monte Carlo simulations, the following conclusions are obtained. First, during a pandemic, the more active the resource support between individuals, the more effectively the disease can be controlled; that is, there is a smaller infection density and a larger epidemic threshold. Second, the disease can be better suppressed when individuals with small degrees are preferentially protected. In addition, there is a critical parameter of contact preference at which the effectiveness of disease control is the worst. Third, the inter-layer degree correlation has a "double-edged sword" effect on spreading dynamics. In other words, when there is a relatively lower infection rate, the epidemic threshold can be raised by increasing the positive correlation. By contrast, the infection density can be reduced by increasing the negative correlation. Finally, the infection density decreases when raising the relative weight of the global information, which indicates that global information about the epidemic state is more efficient for disease control than local information.

A modified Susceptible-Infected-Recovered model for observed under-reported incidence data.

Fitting Susceptible-Infected-Recovered (SIR) models to incidence data is problematic when not all infected individuals are reported. Assuming an underlying SIR model with general but known distribution for the time to recovery, this paper derives the implied differential-integral equations for observed incidence data when a fixed fraction of newly infected individuals are not observed. The parameters of the resulting system of differential equations are identifiable. Using these differential equations, we develop a stochastic model for the conditional distribution of current disease incidence given the entire past history of reported cases. We estimate the model parameters using Bayesian Markov Chain Monte-Carlo sampling of the posterior distribution. We use our model to estimate the transmission rate and fraction of asymptomatic individuals for the current Coronavirus 2019 outbreak in eight American Countries: the United States of America, Brazil, Mexico, Argentina, Chile, Colombia, Peru, and Panama, from January 2020 to May 2021. Our analysis reveals that the fraction of reported cases varies across all countries. For example, the reported incidence fraction for the United States of America varies from 0.3 to 0.6, while for Brazil it varies from 0.2 to 0.4.

Bistability in deterministic and stochastic SLIAR-type models with imperfect and waning vaccine protection.

Various vaccines have been approved for use to combat COVID-19 that offer imperfect immunity and could furthermore wane over time. We analyze the effect of vaccination in an SLIARS model with demography by adding a compartment for vaccinated individuals and considering disease-induced death, imperfect and waning vaccination protection as well as waning infections-acquired immunity. When analyzed as systems of ordinary differential equations, the model is proven to admit a backward bifurcation. A continuous time Markov chain (CTMC) version of the model is simulated numerically and compared to the results of branching process approximations. While the CTMC model detects the presence of the backward bifurcation, the branching process approximation does not. The special case of an SVIRS model is shown to have the same properties.

Bayesian inference for asymptomatic COVID-19 infection rates.

To strengthen inferences meta-analyses are commonly used to summarize information from a set of independent studies. In some cases, though, the data may not satisfy the assumptions underlying the meta-analysis. Using three Bayesian methods that have a more general structure than the common meta-analytic ones, we can show the extent and nature of the pooling that is justified statistically. In this article, we reanalyze data from several reviews whose objective is to make inference about the COVID-19 asymptomatic infection rate. When it is unlikely that all of the true effect sizes come from a single source researchers should be cautious about pooling the data from all of the studies. Our findings and methodology are applicable to other COVID-19 outcome variables, and more generally.

Inhibition of the hexamerization of SARS-CoV-2 endoribonuclease and modeling of RNA structures bound to the hexamer.

Non-structural protein 15 (Nsp15) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) forms a homo hexamer and functions as an endoribonuclease. Here, we propose that Nsp15 activity may be inhibited by preventing its hexamerization through drug binding. We first explored the stable conformation of the Nsp15 monomer as the global free energy minimum conformation in the free energy landscape using a combination of parallel cascade selection molecular dynamics (PaCS-MD) and the Markov state model (MSM), and found that the Nsp15 monomer forms a more open conformation with larger druggable pockets on the surface. Targeting the pockets with high druggability scores, we conducted ligand docking and identified compounds that tightly bind to the Nsp15 monomer. The top poses with Nsp15 were subjected to binding free energy calculations by dissociation PaCS-MD and MSM (dPaCS-MD/MSM), indicating the stability of the complexes. One of the identified pockets, which is distinctively bound by inosine analogues, may be an alternative binding site to stabilize viral RNA binding and/or an alternative catalytic site. We constructed a stable RNA structure model bound to both UTP and alternative binding sites, providing a reasonable proposed model of the Nsp15/RNA complex.

Prediction of dengue fever outbreaks using climate variability and Markov chain Monte Carlo techniques in a stochastic susceptible-infected-removed model.

The recent increase in the global incidence of dengue fever resulted in over 2.7 million cases in Latin America and many cases in Southeast Asia and has warranted the development and application of early warning systems (EWS) for futuristic outbreak prediction. EWS pertaining to dengue outbreaks is imperative; given the fact that dengue is linked to environmental factors owing to its dominance in the tropics. Prediction is an integral part of EWS, which is dependent on several factors, in particular, climate, geography, and environmental factors. In this study, we explore the role of increased susceptibility to a DENV serotype and climate variability in developing novel predictive models by analyzing RT-PCR and DENV-IgM confirmed cases in Singapore and Honduras, which reported high dengue incidence in 2019 and 2020, respectively. A random-sampling-based susceptible-infected-removed (SIR) model was used to obtain estimates of the susceptible fraction for modeling the dengue epidemic, in addition to the Bayesian Markov Chain Monte Carlo (MCMC) technique that was used to fit the model to Singapore and Honduras case report data from 2012 to 2020. Regression techniques were used to implement climate variability in two methods: a climate-based model, based on individual climate variables, and a seasonal model, based on trigonometrically varying transmission rates. The seasonal model accounted for 98.5% and 92.8% of the variance in case count in the 2020 Singapore and 2019 Honduras outbreaks, respectively. The climate model accounted for 75.3% and 68.3% of the variance in Singapore and Honduras outbreaks respectively, besides accounting for 75.4% of the variance in the major 2013 Singapore outbreak, 71.5% of the variance in the 2019 Singapore outbreak, and over 70% of the variance in 2015 and 2016 Honduras outbreaks. The seasonal model accounted for 14.2% and 83.1% of the variance in the 2013 and 2019 Singapore outbreaks, respectively, in addition to 91% and 59.5% of the variance in the 2015 and 2016 Honduras outbreaks, respectively. Autocorrelation lag tests showed that the climate model exhibited better prediction dynamics for Singapore outbreaks during the dry season from May to August and in the rainy season from June to October in Honduras. After incorporation of susceptible fractions, the seasonal model exhibited higher accuracy in predicting outbreaks of higher case magnitude, including those of the 2019-2020 dengue epidemic, in comparison to the climate model, which was more accurate in outbreaks of smaller magnitude. Such modeling studies could be further performed in various outbreaks, such as the ongoing COVID-19 pandemic to understand the outbreak dynamics and predict the occurrence of future outbreaks.

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference.

We present a two-headed approach called Bayesian Integrated Coalescent Epoch PlotS (BICEPS) for efficient inference of coalescent epoch models. Firstly, we integrate out population size parameters, and secondly, we introduce a set of more powerful Markov chain Monte Carlo (MCMC) proposals for flexing and stretching trees. Even though population sizes are integrated out and not explicitly sampled through MCMC, we are still able to generate samples from the population size posteriors. This allows demographic reconstruction through time and estimating the timing and magnitude of population bottlenecks and full population histories. Altogether, BICEPS can be considered a more muscular version of the popular Bayesian skyline model. We demonstrate its power and correctness by a well-calibrated simulation study. Furthermore, we demonstrate with an application to SARS-CoV-2 genomic data that some analyses that have trouble converging with the traditional Bayesian skyline prior and standard MCMC proposals can do well with the BICEPS approach. BICEPS is available as open-source package for BEAST 2 under GPL license and has a user-friendly graphical user interface.[Bayesian phylogenetics; BEAST 2; BICEPS; coalescent model.].

Sensitivity of SARS-CoV-2 Life Cycle to IFN Effects and ACE2 Binding Unveiled with a Stochastic Model.

Mathematical modelling of infection processes in cells is of fundamental interest. It helps to understand the SARS-CoV-2 dynamics in detail and can be useful to define the vulnerability steps targeted by antiviral treatments. We previously developed a deterministic mathematical model of the SARS-CoV-2 life cycle in a single cell. Despite answering many questions, it certainly cannot accurately account for the stochastic nature of an infection process caused by natural fluctuation in reaction kinetics and the small abundance of participating components in a single cell. In the present work, this deterministic model is transformed into a stochastic one based on a Markov Chain Monte Carlo (MCMC) method. This model is employed to compute statistical characteristics of the SARS-CoV-2 life cycle including the probability for a non-degenerate infection process. Varying parameters of the model enables us to unveil the inhibitory effects of IFN and the effects of the ACE2 binding affinity. The simulation results show that the type I IFN response has a very strong effect on inhibition of the total viral progeny whereas the effect of a 10-fold variation of the binding rate to ACE2 turns out to be negligible for the probability of infection and viral production.

The evolution of happiness pre and peri-COVID-19: A Markov Switching Dynamic Regression Model.

Happiness levels often fluctuate from one day to the next, and an exogenous shock such as a pandemic can likely disrupt pre-existing happiness dynamics. This paper fits a Marko Switching Dynamic Regression Model (MSDR) to better understand the dynamic patterns of happiness levels before and during a pandemic. The estimated parameters from the MSDR model include each state's mean and duration, volatility and transition probabilities. Once these parameters have been estimated, we use the one-step method to predict the unobserved states' evolution over time. This gives us unique insights into the evolution of happiness. Furthermore, as maximising happiness is a policy priority, we determine the factors that can contribute to the probability of increasing happiness levels. We empirically test these models using New Zealand's daily happiness data for May 2019 -November 2020. The results show that New Zealand seems to have two regimes, an unhappy and happy regime. In 2019 the happy regime dominated; thus, the probability of being unhappy in the next time period (day) occurred less frequently, whereas the opposite is true for 2020. The higher frequency of time periods with a probability of being unhappy in 2020 mostly correspond to pandemic events. Lastly, we find the factors positively and significantly related to the probability of being happy after lockdown to be jobseeker support payments and international travel. On the other hand, lack of mobility is significantly and negatively related to the probability of being happy.

Improved log-Gaussian approximation for over-dispersed Poisson regression: Application to spatial analysis of COVID-19.

In the era of open data, Poisson and other count regression models are increasingly important. Still, conventional Poisson regression has remaining issues in terms of identifiability and computational efficiency. Especially, due to an identification problem, Poisson regression can be unstable for small samples with many zeros. Provided this, we develop a closed-form inference for an over-dispersed Poisson regression including Poisson additive mixed models. The approach is derived via mode-based log-Gaussian approximation. The resulting method is fast, practical, and free from the identification problem. Monte Carlo experiments demonstrate that the estimation error of the proposed method is a considerably smaller estimation error than the closed-form alternatives and as small as the usual Poisson regressions. For counts with many zeros, our approximation has better estimation accuracy than conventional Poisson regression. We obtained similar results in the case of Poisson additive mixed modeling considering spatial or group effects. The developed method was applied for analyzing COVID-19 data in Japan. This result suggests that influences of pedestrian density, age, and other factors on the number of cases change over periods.

Hotspots for mutations in the SARS-CoV-2 spike glycoprotein: a correspondence analysis.

Spike glycoprotein (Sgp) is liable for binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the host receptors. Since Sgp is the main target for vaccine and drug designing, elucidating its mutation pattern could help in this regard. This study is aimed at investigating the correspondence of specific residues to the SgpSARS-CoV-2 functionality by explorative interpretation of sequence alignments. Centrality analysis of the Sgp dissects the importance of these residues in the interaction network of the RBD-ACE2 (receptor-binding domain) complex and furin cleavage site. Correspondence of RBD to threonine500 and asparagine501 and furin cleavage site to glutamine675, glutamine677, threonine678, and alanine684 was observed; all residues are exactly located at the interaction interfaces. The harmonious location of residues dictates the RBD binding property and the flexibility, hydrophobicity, and accessibility of the furin cleavage site. These species-specific residues can be assumed as real targets of evolution, while other substitutions tend to support them. Moreover, all these residues are parts of experimentally identified epitopes. Therefore, their substitution may affect vaccine efficacy. Higher rate of RBD maintenance than furin cleavage site was predicted. The accumulation of substitutions reinforces the probability of the multi-host circulation of the virus and emphasizes the enduring evolutionary events.